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Abstract 

Non-equilibrium steady states of quantum fields on star graphs arc explicitly con- 
structed. These states are parametrized by the temperature and the chemical potential, 
associated with each edge of the graph. Time reversal invariance is spontaneously bro- 
ken. We study in this general framework the transport properties of the Schrodinger 
and the Dirac systems on a star graph, modeling a quantum wire junction. The inter- 
action, which drives the system away from equilibrium, is localized in the vertex of the 
graph. All point-like vertex interactions, giving rise to self-adjoint Hamiltonians possibly 
involving the minimal coupling to a static electromagnetic field in the ambient space, are 
considered. In this context we compute the exact electric steady current and the non- 
equilibrium charge density. We investigate also the heat transport and derive the Casimir 
energy density away from equilibrium. The appearance of Friedel type oscillations of the 
charge and energy densities along the edges of the graph is established. We focus finally 
on the noise power and discuss the non-trivial impact of the point-like interactions on the 
noise. 
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1 Introduction 



Transport properties of quantum wire junctions attracted in the last two decades much attention 
[T]-[52]. The experimental realizations of quantum wires include nowadays carbon nanotubes, 
semiconductor, metallic and polymer nanowires, and quantum Hall edges. While the equilib- 
rium features of these devises has been extensively explored, there is recently a growing interest 
in non-equilibrium phenomena. A typical problem in this context is schematically represented 
in Fig. [1} A quantum wire junction is modeled by a star graph T with n edges (leads) Ei, 
each of them connected to a heat reservoir (bath) with (inverse) temperature Pi and chemical 
potential /ij. Assume for simplicity that Ei are infinite and that the interaction, described by 
a scattering matrix §(fc), is localized at the vertex V of F. The system is away from equilib- 
rium if Ei{k) admits at least one non-trivial transmission coefficient among edges with different 
temperature and/or chemical potential. Two main questions arise at this point. The first one 
concerns the existence of a steady state describing the above situation. Provided that such a 
state exists, it is natural to ask about the general features of the system in this state. In the 
present paper we show that the first question has an affirmative answer, constructing explicitly 
a suitable steady state flp,^, which is parametrized by /3 = (/3i, ...,/3n) and /i = (/ii, ...,/i2) and 
captures the non-equilibrium properties of the system. Afterwards we derive the expectation 
values of several basic physical observables (currents, charge and energy densities,...) in ^2/3,/^, 
which characterize the system and thus answer the second question. 




Figure 1: A star graph F with scattering matrix §(fc) at the vertex V and edges Ei connected 
at infinity to thermal reservoirs with inverse temperature Pi and chemical potential fii. 

It is perhaps useful to recall that the systems admitting non-equilibrium steady states 
(NESS's) [23] represent an important subclass of the large family of non-equilibrium systems. 
Unlike in equilibrium, a system in a NESS admits steady currents. Nevertheless, all macroscopic 
properties are still time independent like in equilibrium. 

The formulation of a suitable statistical mechanical framework for treating NESS's has not 
been yet completed and is intensively investigated [3l]-|in]. We follow below a microscopic 
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approach and address the question using quantum field theory methods. Our construction of 
NESS generahzes that of a Gibbs state in finite temperature quantum field theory. We develop 
first a general algebraic framework, which applies for any dispersion relation, thus covering both 
the relativistic and non-relativistic cases. As already mentioned, we consider for simplicity fields 
which propagate freely in the bulk F \ of the graph, the interaction being localized exclusively 
in the vertex V. We discuss all possible point-like interactions for which the Hamiltonian of 
our system, being a Hermitian operator in the bulk F \ extends to a self-adjoint operator on 
the whole graph F. We use at this point of our construction some simple results |ll]-|l3] from 
the spectral theory of differential operators on graphs (known also as "quantum graphs"). 

It is worth mentioning that non-equilibrium systems of the type shown in Fig. [T] have been 
investigated in the past by two different methods. The first one is the scattering approach 
initiated by Landauer [H] and completed by Biittiker [15]. The second one is based on the 
linear response theory [H] , originally developed by Kubo for macroscopic samples and extended 
later by Baranger and Stone [H] to mesoscopic systems. The equivalence between these two 
approaches has been demonstrated in Ref. [IS]. 

The construction of the NESS developed in this paper, is formulated in purely alge- 
braic terms and adopts a deformation of the algebra of canonical (anti) commutation relations. 
No approximations, like linear response theory, are involved. The states allow for a unified 
description of systems with different dynamics. We illustrate this fact by treating below both 
the Schrodinger and Dirac equations in the same way. The abstract construction of fi^^^ is 
first tested by reproducing the Landauer-Biittiker (L-B) steady current. Afterwards, taking 
advantage of the exact solvability of the point-like interactions in the state we investigate 
several new features of the non-equilibrium dynamics, applying the general framework to other 
physical observables. We compute the charge density in fi^^^ and show the presence of char- 
acteristic Friedel type oscillations along the leads. We also present a fully microscopic-based 
calculation of the Casimir effect in a quantum wire junction away from equilibrium. In this 
context we establish the non-equilibrium Stefan-Boltzmann law and derive the heat (energy) 
transport. We also compute the exact two-point current-current correlation function away from 
equilibrium and investigate the noise power. In particular, we discuss the impact of the point- 
like interaction in the junction on the behavior of both thermal and shot noises. Finally, we 
generalize all results in the presence of a static electromagnetic field in the ambient space. 

The paper is organized as follows. In the next section we develop a simple algebraic frame- 
work for the the construction of steady states of quantum wire junctions. Section 3 is devoted to 
the Schrodinger junction. In section 4 we change the dynamics, considering the Dirac equation 
and the relative steady states on the junction. The novelty here is the presence of antiparticles 
which contribution and role are considered in detail. In section 5 we discuss the results and 
comment on the possible further developments. 



2 Algebraic construction of the states flp^^ 

Previous investigations [IE]- [21], [IH] have shown that a convenient coordinate system for ex- 
pressing quantum fields on star graphs is provided by the deformed algebras A± of canonical 
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(anti) commutation relations, generated by {ai{k), a*{k) : k G 



, n} which satisfy 



[ai{k) , aj{p)]± = [a*{k) , a* 
[a,{k) , a*{p)]± = 27r[6{k - p)6,, + S^,{k)6{k + p)] 



and the constraints 



,ik) = J2%ik)a,i-k) 



(2.1) 
(2.2) 

(2.3) 



The index ± in (2.1 , 2.2) refers to Fermi/Bose statistics and §(fc) is the nxn scattering matrix 
describing the point-hke interaction in the vertex of the graph. We assume in what follows 
unitarity 

§{k)§{ky = i (2.4) 

and Hermitian analyticity [50]-[53] 

S{k)* = S{-k). (2.5) 

The latter implies that the ^-operation is a conjugation in A±. Combining (2.4) and (2.5) one 
concludes that S(fc)S(— /c) = I, ensuring the consistency of the constraints (2.3). 

Our main goal now is to construct state, (continuous) linear func- 

tional over A±. We recall in this respect that A± is a simplified version of the so called 
reflection-transmission algebra [50]- [53], describing factorized scattering in integrable models 



with point-like defects in one dimension. The Fock and the (grand canonical) Gibbs state over 
A± describe equilibrium physics and have been largely explored [50j-i53j. The physical input 
for constructing the new states fi^^^ is the observation that the sub-algebras A±. and A'^*', 
generated by {ai{k), a*{k) : < 0} and {ai(/c), a*{k) : > 0} respectively, parametrize the 
asymptotic incoming and outgoing fields. Accordingly, both ^™ and A^^ are conventional 
canonical (anti) commutation relation algebras. Indeed, the 5{k + p) term in (2.2) vanishes if 
both momenta are negative or positive. Notice also that the constraints (2.3) relate A^l with 



A°^^. It is crucial for what follows that the whole reflection-transmission algebra A± can be gen- 



erated either by A. 

on this kind of asymptotic completeness property. In fact, we will start with an equilibrium 



or by A^ via (2.3). Our strategy for constructing the NESS ^p^n is based 



state on A^^ and extend it by means of (2.3) to a non-equilibrium state on A±. 



The first step is to describe the asymptotic dynamics and symmetries at t = — oo (i.e. before 
the interaction) in terms of A™. Since the asymptotic fields are free, it is natural to introduce 
the edge Hamiltonians 

dk 

hr= I —uj^{k)a*{k)a^{k) (2.6) 



27r 



and edge charges 



dA; 
2^ 



a*{k)ai{k) 



where uJiik) > is the dispersion relation in the edge Ei. At this point we define 

n 



(2.7) 



(2.8) 
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and introduce the equilibrium Gibbs state over A'^ in the standard way [54j. For any polynomial 
V over A±. we set 

in^,,,V{a:{k,),a,{p,))Q^,,) ^ {V{a:{k),a,{p,)))p,, = ^Tr [e-''V{a*{k,),a,{p,))] , (2.9) 

where ki < 0, pj < and Z = Tr (e~^) . It is well known [54J that one can compute the 
expectation values (2.9) by purely algebraic manipulations and that all these expectation values 
can be expressed in terms of the two-point functions, which are written in terms of the familiar 
Fermi/Bose distributions in the following way 



{aVp)ai{k))p, 



{ai{k)a*Ap))p 



^-I3i[u)i{k)-Hi] 
1 

~ I ± Q-l3iliOi{k)-fii] 



75ij27T5{k -p) , 



Sij27i5{k — p) . 



(2.10) 
(2.11) 



We stress that ( 2.10|2.11 ) hold on i.e. only for negative momenta. 

The second step is to extend (2.9 2.11 ) to the whole algebra A±, namely to positive momenta. 
For this purpose we use the relations (2.3). One finds in this way 

n 

{a*{p)a,{k))p,, = 27r| [di-k)df{k)5,, + d{k) Su{k) df{-k) §,,(-fc)] 5{k - p) 



1=1 



+ 



e{-k)df{k)%{k) + e{k)^,.j{k)df{-k) 6{k + p) 



were for simplifying the notation we introduced 



df{k) 



\ -|- Q- Pi[LOi(k)- pLi] 



The explicit expression of {ai{k)a*{p)) i^^^ is obtained from (2.12) by the substitution 

df{k)^cf{k) 



1 



(2.12) 



(2.13) 



(2.14) 



I ^ Q-Pi[uJi{k)-fj.i] 

The final step is to compute a generic correlation function. Employing the commutation 
relations ( 2.1|2.2 ), one can reduce the problem to the evaluation of correlators of the form 

M N 

n* 



m=l n=l 

which can be computed in turn by iteration via 

M N M 

(11 ('irn{k^JY[a*^"{PjJ)p,^. = Smn J2^ai,{k,,)a*^"^{pjJ)^^^{Y[ai^{k,J JJ a*^"(p,J);3,^ . 

(2.16) 



m=l 



n=l 



m=l 



M 



m=2 



N 



n=l 
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In conclusion, we emphasize once more that the use of the deformed algebras A± in the 
construction of represents only a convenient choice of coordinates, which has the advantage 
to be universal and to apply to a variety of systems characterized by a scattering matrix S(/c). 
In support of this statement, we consider below the Schrodinger and the Dirac equations on 
the star graph F. 



3 The Schrodinger junction 
3.1 Preliminaries 



In this section we apply the general algebraic construction of the state to the Schrodinger 
system on a star graph F with point-like interactions in the vertex V of F. We will consider 
Fermi statistics and for simplifying the notation will omit the apex + in the Dirac distributions 
df{k) and cf{k). As observed in section 3.5 below, most of the results can be easily extended 
to Bose statistics. 

We start by summarizing the main features [18] of the Schrodinger equation on F, recalling 
the description of all point-like interactions leading to a self-adjoint Hamiltonian. Each point 
P in the bulk of F is parametrized by (x, i), where a; > is the distance of P from the vertex V 
and i labels the edge. In the bulk F \ of the graph the Schrodinger field ilj{t, x, i) with Fermi 
statistics satisfies 

idt + ^di^i^{t,x,i) = Q, (3.17) 
with standard equal-time canonical anticommutation relations 

[^(0,Xi,^i) , V'(0,X2,«2)]+ = [V^*(0,Xi,2i) , V^*(0,X2,22)]+ = 0, (3.18) 
[^/'(O, Xi, ii) , V^*(0, X2, ^2)]+ = ^iiia ^{xi - X2) ■ (3.19) 



The interaction in the vertex is fixed by requiring that the bulk Hamiltonian defined by (3.17) 
(essentially the operator —(9^) admits a self-adjoint extension on the whole graph. According 
to some elementary results from the spectral theory |1T]-P3] of differential operators on graphs, 
all such interactions are described by the boundary conditions 

n 

lim V [A(I - - i(I + V)ijd,] ij{t, x,j) = , (3.20) 

x— >0+ •^■"^ 

where U is an arbitrary n x n unitary matrix and A G M is a parameter with dimension of 



mass. Eq. (3.20) guaranties unitary time evolution of the system and generalizes to the graph 
F the familiar mixed (Robin) boundary condition on the half- line M+. The matrices U = I and 
U = — I define the Neumann and Dirichlet boundary conditions respectively. 

The explicit form of the scattering matrix, expressed in terms of U and A, is ||4T]-|^ 

[A(I - U) + MI + U)] ^"^-^^^ 
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and has a simple physical interpretation: the diagonal element Sii{k) represents the reflection 
amplitude from the vertex on the edge Ei, whereas with i ^ j equals the transmission 

amplitude from to Ej. One easily verifies that (3.21) satisfies (2.4, 2.5) and therefore defines 
an algebra A+ of the type introduced in the previous section. Notice that 

§(A)=U, S{-X)=V-\ (3.22) 



showing that the unitary matrix U entering the boundary conditions (3.20) is actually the 
scattering matrix at scale A. 

A remarkable property of (3.21) is that it can be diagonalized for any A; by a /c-independent 

(3.23) 



unitary matrix. In fact, let U be the unitary matrix diagonalizing U, namely 



W-i = = diag (e^^^S e2^°2 



) 



By means of (3.21) one concludes that U diagonalizes Ei{k) for any k as well, and that 

'k + irji k + ir]2 k + i?7„ ~ 



where 



Sd{k) = U*S{k)U = diag 



Tji = Atan(ai 



k -irji k- 1772' k - irj^ 



TT 71 

- < ai< - 
2 - 2 



(3.24) 



(3.25) 



Illustrating various aspects of the Schrodinger junction, we will use below the most general 
2x2 S-matrix 



S{k) 



k^+ik(rii—ri2) cos(g)+r?ir?2 
{fc-i»7i)(A;-iJ72) 

-ic~"^fc(r?i->72) sin(6') k'^—ik{r}i—ri2) cos(g)+qir?2 
(k-irii){k-iri2) {k-irii)(k-iri2) 



— io"^fc(r;i— r;2) sin(0) 
{fc-i»7i)(A:-r^2) 



(3.26) 



where ip and 6 are arbitrary angles. 

The general representation (3.24) implies that S{k) is a meromorphic function in the complex 
/c-plane with finite number of simple poles on the imaginary axis. For simplicity we consider 
in this paper the case without bound states (poles in the upper half plane), referring for the 
general case to [33], |2H], [SI] and the comments in section 3.5. In other words, we assume that 



dk 
2^ 



Akx 



Sij{k) = 0, x>0. 



(3.27) 



It has been shown in Ref. [IBJ that in this case the solution of equation (|3.17|) is fixed uniquely 

)rm 

ai(A;)e-'"(^')*+''^", 



by (3.18 3.20) and takes the following simple form 

""^ dk 



271 



u{k) 



2m 



The equation of motion (3.17) is invariant under the time reversal operation 

Tip{t, X, i)T'^ = -r]Ti'{-t, x,i), \r]T\ = I , 



(3.28) 



(3.29) 
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T being an anti-unitary operator. We stress however that the boundary condition (3.20) pre- 
serves the time reversal symmetry only if U is symmetric [30] , namely 

U* = U. 



The electric current 



jx{t,x,i) = i— [ip*{d:r'ip) - [dx'ipl'ip] {t,x,i) 



(3.30) 
(3.31) 



e being the electric charge, and the energy current 

ext{t,x,t) = ^[(dtr) (d.^P) + (d.r) (dt^) - (dtd^r)^- r (9A^)](t,x,0 , (3.32) 
are among the basic physical observables. The time components of these currents are 

jt{t,x,i) = e{ip*'ilj){t,x,i) , (3.33) 
1 



respectively and, as a consequence of (3.17), satisfy the local conservation laws 

{dtjt - dxjx) (t, X, i) = {dtdtt - d^d^t) {t, x,i) = 0. 



Equations (3.35), combined with the Kirchhoff rules 

n n 

J]j,(t,0,z) = 0, J]^^,,(t,0,z) = 0. 



(3.34) 
(3.35) 

(3.36) 



1=1 



ensure the charge and energy conservation in the system. Since the proof of (3.36) at the 
quantum level is a quite subtle, we provide the main steps, focussing for instance on the electric 
current. Using the basic definitions, one easily derives the representation 

-.[w{k)-u{p)] 



X 



e dk dp =,r 

' ' ^ 2m 27r J^^ 27r 

Yl 0) {p; 0) - [^xX*^] {k; 0)xa{p; 0)]ai{p) 



(3.37) 



where 



X{k; x) = e'*^"I + e"^^" §* (k) , x* {k; x) = e"''^"! + e'^" §(fc) . (3.38) 
The trick now is represent the right hand side of ( |3.38[ ) as a boundary term of an integral over 

the half line, namely 



° dk r° dp 



2m / ^ 2tt / ^ 2tt 



oo ^'^ Jo 



Ya*j{k)\x*jiik;x) [dlxu] {p;x) - [dlx*ji] {k; x)xuip; x)\ai{p) 



X 

^ £ '^^'^''^'^-^^'Hk' - P') dxj^^^^ (3.39) 



The final step is to apply the orthogonahty relations 



i=i 



(3.40) 



which hold [E] for § given by (3.21). 

At this stage we are ready to investigate the properties of the system in the state Being 



expressed in terms of the algebra A+, the solution (3.28) and the observables (3.31 3.34) apply 



for any representation of this algebra. This fundamental universality property has been already 
largely explored in the Fock and the Gibbs representations of ^4., which describe equilibrium 
physics. In order to study the non-equilibrium properties of the Schrodinger system in Fig. [T| 
we apply below the representation generated by the state f^/s,^. Since antiparticle excitations 
are absent, we assume without loss of generality that /ii > 0. 



There are two non-trivial two-point correlation functions. Using (2.12), one finds 



' Jo 2n 
dk 



\kx12 



1=1 



(3.41) 



(V^(ti,xi,i)V^*(t2,a;2, j))/3. 



-iw(k)t\2 



27r 



5,,Q(fc)e^'=^- + §i,-(A;)c,(fc)e''=^- + %{k)c,{k)e-'''~^'^ + ^ El,{k)ci{k)%{k)e-'^^^^ 



1=1 



(3.42) 



where ti2 = ti — ^2; 2:12 = Xi — X2 and X12 = Xi+ X2- The invariance of (3.41 , 3.42) under time 



translation implies energy conservation. For systems away from equilibrium one expects that 



the time reversal (3.29) symmetry is instead broken. In fact. 



{'4)*{ti,xi,i)i){t2,X2,j))p,t, ^ (V^*(-t2,a;2, j)V^(-ti,xi,i))/3,^. (3.43) 
showing that TQ^^^ ^ ^/s,^, i-e. that fi^j^^ breaks down spontaneously time reversal invariance 



even if (3.30) holds 



3.2 Charge transport and density 



Using the explicit form of the correlation function ( |3.41 ) and applying standard point-splitting 
technique, one obtains 

Ji(/3,/i) = {i^{t,x,i))p^^ = 

g 

^2m KV^*(^i>a;i,i)<9x2V^(t2,a;2,0)/3,M - {dc,,ip*{ti,xi,i)ilj{t2,X2,i))i3,,,] = 



e I ™ dk 



m In 2tt 



kJ2[Sij-Kik)\']dAk) 



(3.44) 



where is given by (3.21) and covers all point-like interactions leading to a self-adjoint 
Schrodinger Hamiltonian on the star graph. The current Ji{(3,fi) is t- independent, implying 
that are indeed steady states. We emphasize that (3.44) is precisely the Biittiker multi- 
channel generalization [IS] of the Landauer expression [44J for the steady current. It is also 



worth stressing that in our context (3.44) is an exact formula, which does not relay on linear 
response theory. These remarkable features of the states fi^^^ confirm their physical relevance 
and suggest to call them L-B states. 

Let us summarize some of the basic properties of the steady current Jj(/3,/i). First of all 



Jj(/3,/i) is homogeneous (x-independent) and, because of unitarity (2.4), satisfies 



5^J,(/5,/i) = 



(3.45) 



i=l 



This is the manifestation of the operator Kirchhoff rule (3.36) and represents a non-trivial 



check on (|3.44|). Moreover, there are two particular cases in which the system is in equilibrium 

must therefore vanish. The first one is when all thermal reservoirs are 
■ ■ = (3n and /Ui = /i2 = ■ ■ ■ = /in)- In fact, (2.4) implies in this case 

Ji(/3,/i) = 0. (3.46) 



/32 



Another possibility to be in equilibrium is when all transmission coefficients vanish and the 
leads are therefore isolated. In this case §(A;) is diagonal. 



(3.47) 



which implies (3.46) as well. 



In order to illustrate the role of the §-matrices (3.21), it is instructive to consider (3.44) for 
n = 2. Using the general expression (3.26), one finds 



Ji(/3,yu) = -J2(/3,/x) 



m 



■livi 



r]2 ) sm 



dA; 



27r(A;2 + r/2)(A;2 + ^2) 



[di{k) - d2ik)] , (3.4J 



where the sign-difference between Ji and J2 reflects the orientation of the leads. The k- 



integration in (3.48) can not be performed in a closed analytic form, but the integral is well- 



defined and can be computed numerically. The contour plots of Ji for fixed e, 9 and m, displayed 



in Fig. 2, give an idea about the behavior of (3.48 ) in the variables (/3i, /32), (/^i, /i2) and (r^i, 772). 
As usual, higher regions are shown in lighter shades. The plot on the left concerns Ji in the 
plane (32) for fixed {fii, ^2) and (7/1,772). The plot in the middle illustrates the behavior of 
Ji as a function of (/ii,/X2), the variables (/3i,/32) and (771,7/2) being fixed. Finally, the plot on 
the right shows the dependence on the §-matrix variables (7/1,7/2) at fixed temperatures and 
chemical potentials. 

The plots in Fig. 3 are obtained from those in Fig. 2 by fixing /3i, /ii and 7/1 respectively. 
The sign change of Ji in Fig. 3 indicates that varying (32, /i2 and 7/2 one can invert the direction 
of the current. 
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Figure 2: Contour plots of the current Ji in the plane (/3i,/32), (/ii,/i2) and (771,772) respectively. 






Figure 3: Plots of Ji as a function of /32, /i2 and 772 respectively, with all other variables fixed. 



The expression (3.44) significantly simplifies at criticality, i.e. for scale invariant point-like 

(3.49) 



interactions in the vertex of F. In this case the scattering matrix takes the form 
9 being the Heaviside step function. Plugging ( |3.49[ ) in (3.44), one gets 



(3.50) 



i=i 



The high and zero temperature limits of (3.50) are 



Ji(0,/i) = lim Ji{/3,fi) = ^Y1 i^^^ ~ 



Ji(oo,/i) = lim Ji(/3,/i) = 7^Y] - |Uijf ) /ij 



respectively and are related by 



(3.51) 



(3.52) 



(3.53) 
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By means of (3.44) one can derive the conductance tensor Gjj(/3,/i) defined by 



J,(/3,/x) = 5^G,,(/3,/i)V^, 



] 1 



(3.54) 



wliere Vj is the voltage apphed at the edge Ej. It is well known ^56j that Gij{f3,fi) depends 
on the point on Ej where the voltage is applied. Assuming that this point is deeply in the 
reservoir with chemical potential Hj, one has [56] 



Combining (3.44), (3.54) and (3.55) one obtains 

-2 roo ^ 



,(/3,/x) 



m Jq 2ti jjj 



which satisfies Kirchhoff 's rule 



^G,,(/3,/i) = 



(3.55) 



(3.56) 



(3.57) 



i=l 



as it should be. 

Let us focus now on the charge density distribution {jt{t, x, i))i3,^ in the state ^/s,^, which can 
be computed following the above procedure as well. As expected from the current conservation 
(3.35), the charge density 



°° dA; 
2^ 



n 

[§,,(fc)e-2'*^^ + Su{k)e'''^ + 1] d,{k) + ^ \S,,{k)\^d,{k) 



(3.58) 



is time independent as well. There are however two essential novelties with respect to the 
current: 

(a) Qi does not vanish at equilibrium; 

(b) Qi depends on the position x. 

Concerning point (a), we observe that at equilibrium (n isolated leads with §(/c) defined by 
(3.47)) one has 



^•'i(/3,/i,x) = e 



°° dk 
2^ 



{ [S,,{k)e-'"'^ + S,,{k)e^"'^' + 2] d,{k)} . (3.59) 
Therefore, the non-equilibrium charge distribution is x-independent and is given by 

^ z2 l-^'^ ~ dj{k) . (3.60) 
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The x-dependence, mentioned in point (b), is carried by 



f?r(/3,/i;x) 



271 



(3.61) 



which oscillates with the distance x from the vertex. The appearance of such Friedel-type 
oscillations |55j confirms once more that the junction behaves indeed as a point-like defect. 



Since the integration in (3.61) can not be performed in closed form, in order to get an idea 

mperature limit 

[§ii(fc)e-2i'=- + §,,(fc)e2'^^] . (3.62) 



about the oscillations, it is useful to consider the zero-temperature limit 



g°'"' {oo, iJ,;x) = lim g°'"' {jS , /j,; x) = e 

/3fe=/3-s>oo jQ 



27r 



At criticality (3.49) and setting Vu = Va for simplicity, one finds 

eUii 



^r(oo,/i;x) 



sm 



nx 



{2x^/2mfii 



(3.63) 



which shows that the amplitude of the oscillations on the graph decays with the distance from 
the vertex like x~^, which is a typical behavior in one space dimension. 



3.3 Casimir effect away from equilibrium and heat flow 

The structure of the energy density in the state fi/?^^ resembles very much that of the charge 



density (3.58). One has 



£i{x; 13, fi) = {9u(t,x,i)),3,i 



dk 

— c 
2tx 



-UJ 



(k) \ [Su{k)e-^"'^ + Mk)e'''^ + l] d,{k) + \^^Jm%ik) 



(3.64) 



This result confirms the presence of Friedel oscillations in the energy density as well. It is 

ensity 

S,,{k)e^"'^ + 2] d,{k), 



instructive to compare (3.64) to the equilibrium energy density 

dA; 



S-'^{x;/3,fi) 



27r 



toik) [§,,(A;)e 



-2ikx 



associated to (3.47)). One finds 



1 r°° dk 

Sr{x;^,^^)-S,{x;^,^^) = - / -co{k)J2[S^J ' \%m'] dji^) 



2 Jo 27r 



(3.65) 



(3.66) 



which gives the genuine non-equilibrium part of the energy density (3.64). The a;- independent 

(3.67) 



contribution to (3.64), namely 

J 1 •-'0 
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is the Stefan-Boltzmann law in the present context. At criticahty 

where Li^ is the poly logarithm function. 

The counterpart of the L-B formula for the heat (energy) flow is 



(3.68) 



dk 

m In 2tt 



ku{k) J2 - d,{k) = 7-(/3,/i) . (3.69) 



Apart from the additional u{k) factor in the integrand of (3.69), the charge ( |3.44[ ) and energy 
(3.69) flows have the same structure. For this reason 7I(/3,/i) shares with Jj(/3,/i) the general 
properties listed after equation (3.44). In the scale invariant case the energy flow is 

dk fc3g-/3j[w(fc)-A(j] 



7I(/3,/x) 



2m2 



27r 1 + Q-^j['-'{k)-,ij] 



1 1 



(3.70) 



which gives in the zero temperature limit 

1 " 

7i(oo, /i) ^ =-y: (5., - . 



(3.71) 



3.4 Noise 

In this section we derive the noise power generated by the point like interactions in the 
Schrodinger junction. For this purpose we need [57j-[59] the two-point connectec? current-current 
correlator. After some algebra one finds 



{jx{ti,Xi,i)j^{t2,X2,j))f^^'^ = {jx{ti,Xl,i)jx{t2,X2,j))f},f, - {jxih, Xi, i)) i3^^{j^{t2, X2, j)) (3. 



e 



2 



4m2 

>^{xh{ki;xi) [d^Xim] {k2;xi) - [d^Xu] {h; Xi)xim{k2; Xi] 
X {x*mj{k2; X2) [d^Xji] {ki;x2) - [d:,x*mj] ih] X2)xji{ki; ^2)^ , (3.72) 
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where x ^i-nd x* ^"^^ given by (3.38). Using the time translation invariance of (3.72), the noise 
power is defined [59] by 



(3.73) 



P,,(/3,/i;xi,X2;u;) = / dt e'^* 0,.(t, Xi, ^)j,(0, X2, j))™^ 

J — oo 

The zero-frequency hmit {zero-frequency noise power) 

turns out to be Xi^2-independent and is given by: 

,2 roc 



(3.74) 



-k 



5ijdi{k)ci{k) - \Si,{k)\'dj{k)cj{k) - \S,i{k)\'di{k)ci{k) + 



m Jq 2-k 

+ ^ 'Bii{k)ci{k)Eji{k)SjrrXk)drn{k)Eirn{k) 



Lm=l 



(3.75) 



It is instructive to summarize at this point the general features of (3.75): 



(i) Pij{/3,n) is symmetric in i and j. The first three terms of the integrand are manifestly 
symmetric. Concerning the last term, using the identity Ci{k) = 1 — di{k) one gets 



^ Su{k)ci{k)E>ji{k)E>jmik)d,rnik)Simik) = 

l,m=l 

n n 
Sij'^Sim{k)dmik)Eimik) T ^ Sil{k)dl{k)Ejl{k)Sjmik)dm{k)Eim{k) 



(3.76) 



m=l 



Lrn=l 



which is symmetric as well. One can therefore rewrite Pij{P,fi) in the following manifestly 
symmetric form: 



t2 /-oo 



dk. 



Pij{P,fi) = - I ^k\ 6i,d,{k)c,{k) - \Sij{k)\^dj{k)cj{k) - \S,i{k)\^di{k)ci{k) + 

^ Jo ■^'^ 
1 " 

+ - J2 Mk)Mk)^jUk)^iUk)[ci{k)dUk) + cUk)di{k)]\ . (3.77) 



Lm=l 



(ii) The last identity implies also that Pij{P,fi) is real; 

(iii) As expected, Pij{f3,fi) satisfies the Kirchhoff rule 

n n 

J]P,,(/3,/i) = J]P,,(/3,/i) = 

i=i j=i 



(3.78) 
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which provides an useful check. One has actually 



^ Pij{l3, yu; 0, X2; u) = ^ Pij{l3, /i; a;i, 0; w) = , 



(3.79) 



i=l 



at any frequency u. 



(iv) All noise components Pij{l3,jj,) vanish for isolated leads (3.47). 



Let us discuss now the behavior of the noise, starting with the case n = 2. Combining (3.26) 



with (3.77), we find 



Pii(/3,/x) = [e(?7i - r]2) sin 
di{k) + d2ik) ~ 2di{k)d2ik) - 



m 



dk 







27r (P + r72)(P + r72; 



1 -?72jsm 



{k^ + rf^){k^ + ri) 



■[di{k) - d2{k)f 



(3.80) 



For r/i = ri2 and/or = the leads are isolated (see (3.26)) and the noise vanishes according 



to point (iv) above. Like for the steady current, we report some contour plots, showing the 
comphcated dependence of the noise on the parameters (/3i,/32), (yUi,/i2) and (?7i,?72) for fixed 
e, 6' and m. Fig. 4 illustrates the behavior of Pii{P,fi) in each pair of these variables, the 
remaining two being fixed. The left plot is the noise in the plane {/3i,P2)- In the middle we 
display the noise as a function of the chemical potentials (/ii,/X2)- Finally, the right plot shows 
the dependence on the S-matrix variables {r]i,r)2). 






Figure 4: Contour plots of the noise Pu in the plane (/3i,/32), (/Wi,/i2) and {rii,ri2) respectively. 

The plots in Fig. 5, obtained from those in Fig. 4 by fixing Pi, Hi and r)i respectively, confirm 
that the noise Pij{P, /i) depends in a complicated way on /3, fi and the §-matrix parameters. At 
criticality however, as expected on general grounds, the situation simplifies and one can push 
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Figure 5: Plots of Pn as a function of /32, and r/2 respectively, with all other variables fixed. 



further the analytic computation. In fact, inserting (3.49) in (3.77) one gets. 



1 " — — I 



Lm=l 



where 



dk 
2^ 



k di{k)cj{k) 



(3.81^ 



(3.82) 



For equal temperatures Pi = (3j = (3 the integration in (3.82) can be performed exphcitly and 
one finds 

( m cPf^ -r _ — 



Therefore, in the case /3i = 
obtains 



i'^,(/3,/x) 



27r/3 



Pn = P with generic chemical potentials /ij > one 



-|u. 



1 + e/3Mi ' 1 + e/^/^j- 



-)- e/^/^i 



+ 



1 — — e^^' + e'^'''" / 1 + e^^' 



Z ,m — 1 



(3.84) 



It is instructive to consider at this stage the two limits leading to the shot and thermal 

;h are 

(3.85) 



noise. For deriving the shot noise, we need the /3 — )■ oo limit of the integrals (3.83), which are 

if /ij < fij , 

:(/ii - fij) if /ij > yUj . 



hm h,{P) 



m I 



Therefore, in the scale invariant case (3.49) the shot noise is 

,2 " 



PijijJ,) = lim PijiP^jJ,) = UilUjlUjm^inillJ'l -/iri 

/3— !>oo 47r 



(3.86) 



Z ,m— 1 
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which exhibits the standard behavior [57j-[59| in terms of — fij\- 

In order to compute the thermal noise, we consider (3.84) for /ii = fi2 

has 



e2 



In the purely thermal case (/i — )■ 0) one finds 



2nf3 



Hn = /i. One 
(3.87) 

(3.88) 



which is the well-known Johnson-Nyquist formula. 

A remarkable feature of the thermal noise is that away from criticality the point-like inter- 
actions at the vertex can modify the linear behavior (3.88) for large T. Let us consider indeed 
(3.80) for Pi = ^2 = P and /ii = /i2 = 0, namely 

""^ dfc 



Pii(/3) = 2[e(r/i-r72)sin 



271 



(k^ + r]l){k^ + r]i) 1 + e 



fc2 



(3.89) 



Assuming for simplicity 772 = and introducing the variables rj = rji and ^ = e" 



'P2^^ one gets 



^11 (/3) = [er^sin 



1 



2nm 



d^ 



1 



£ - HO 



(3.90) 



Since the fc- integration in (|3.26|) can not be performed exactly, in order to estimate the tem- 



perature dependence one can use the inequalities 



[erj sin(^^)]' 



1 



ATcm \ 2m 



< PiiiP) < [er]sm 



27rm \ 2m 



where 



/ d^l — — = -e'^Ei(-a) 

Jo [a - HOI 

Ei being the exponential integral function. In this way one finds 



1 



\in(r) 



for T ^ , 
for T 00 , 



(3.91) 



(3.92) 



(3.93) 



which shows that the fc-dependence of the §-matrix indeed modifies the Johnson-Nyquist behav- 
ior at high temperatures. The milder logarithmic divergence for large T provides an attractive 
experimental signature. 
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3.5 External electromagnetic field 



In the above considerations the interaction was locahzed in the junction. We extend here the 
framework to the more reahstic physical situation of a junction in a three-dimensional ambient 
space with a classical static magnetic field, interacting with the Schrodinger excitations along 
the leads. The graph F, modeling the junction, is embedded in M^, equipped with a Cartesian 
coordinate system whose origin O coincides with the vertex V oi T . The direction of each 
edge Ei C is determined by the unit vector e^*-*. At any point P G the magnetic field 
B(P) = rot[A(P)] is generated by the potential A(P). The minimal coupling of the Schrodinger 
field if with A gives the following equation of motion 



id. 



1 

2m 



{\d^ - eA^{x, i)) (id^ - eA^{x, i)) 



0, 



Ax{x,i) being the projection 

A,(x,z) = e«-A(P), 



P = (x,i) er c 



(3.94) 



(3.95) 



of the potential A along the edge Ei. All self-adjoint extensions of the relative Hamiltonian are 
now parametrized by the boundary conditions [61] 



lim V [A(I - U),, - (I + lJ)i,{id, - A,{x,j))] (^(t, X, j) = , 



a; -5-0 



(3.96) 



The conserved electric current is 

. e 



j^{t,x,i) = i— [(p*{dx(p) - {d^(p*)(p] {t,x,i) A^(a;, ii)(v3V)(t, x, z) 

Im m 



(3.97) 



It is easy to show now that the solution of the problem ( 3.94[3.96 ) can be reduced to that 
described in section 3.1. Indeed, let us introduce 



ijj(t,x,i)=e "'°'''''''^(p{t,x,i) , a{x,i) 



dyAy(y,i) 



(3.98) 



where we assumed that Ax{x,i) are integrable on the half line. Notice that ip and cp have the 
same behavior for x — )■ oo. Moreover, ip satisfies (3.17) and in terms of ip the current (3.97) 
takes precisely the form (3.31). The interaction is totally absorbed in the boundary condition 
for ijj following from ( 3.96[3.98 ). One has 



lim V [A(I - U(A)),, - i(I + U(A)),,9,] ij{t, x,j) = 

■ i.n+ ' * 



(3.99) 



where 



U.,(A) 



a(0, i) . 



(3.100) 
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Combining (3.21) and (3.100) one concludes that the substitution 



{3.101] 



extends all the results of this section to the case of a junction minimally coupled to a time- 
independent ambient magnetic field. We stress that the correlation functions of fields localized 
in different edges (see e.g. ( 3.41|3.42|3.72 )) are sensitive to the transformation (3.101). In 
particular, the field A has a non-trivial impact on the noise power at frequency u ^ 0. We will 
analyze this issue in more details elsewhere. 



3.6 Remarks 



Let us discuss first the role of possible bound states of which have been excluded in 



the above considerations by assuming (3.27). It has been shown in previous work 



) that the bound states generate new quantum degrees of freedom, which have a non-trivial 



contribution to the correlation function (3.41). The key point is that this contribution depends 



on the space-time coordinates only through the combinations ti2 and Xi2- According to (3.31 



3.32), the charge (3.44) and energy (3.70) flows are therefore not affected by the presence of 



bound states. The relative densities however get [5^ nontrivial bound state contributions. 



One can investigate along the above lines also the Schrodinger equation (3.17) with Bose 
statistics. The final results in this case obviously follow from equations ( 3.41[3.44[3.58 3.70[3.77 ) 
by substituting the Fermi distribution with the Bose distribution d^{k). Concerning the noise, 
in the scale invariant case one gets for bosons 



271(3 



+ 



1 = 1 l.m = l 



^Pi^'i _|_ g/^Mm /I — e^^' 
In 



1 - e^/^' 



(3.102) 



which, compared to (3.84), shows haw the zero frequency noise power depends on the statistics. 



For instance, in the shot noise limit one obtains 



Pj - (/i) = lim Pj . (/3,/i) = UiiUjiUi.mUimlyU; - /^m| , 

p^oo ■' 47r 

i ,m — 1 



(3.103) 



which has the magnitude of the fermionic shot noise (3.86) but the opposite sign |59 
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4 The Dirac junction 



4.1 Preliminaries 

The massless Dirac equation on the star graph F is 

{itdt - lxdx)ip{t, X, i) = . 

where 



i){t,x,i) 



ip2it,x,i) 



It 



1 

1 



X > 



1 
-1 



(4.104) 



(4.105) 



We assume that ipa satisfy the conventional equal-time anti-commutation relations. The bound- 
ary conditions which define all self-adjoint extensions of the bulk Hamiltonian i'^tlxdx are |62j . 



(4.106) 



^i(t,o,o = 5];%^2(t,o,j), 

where U is any unitary n x n matrix. In physical terms U parametrizes all point-like interactions 
for which ijt'Jxdx extends to a self-adjoint Hamiltonian to the whole F. Observing that both 



the equation of motion (4.104) and the boundary condition (4.106) preserve scale invariance, it 



is not surprising that the scattering matrix corresponding to these interactions is simply (3.49). 
The Dirac field ip is complex, has a relativistic dispersion relation 



u{k) = \k\ 



(4.107) 



and describes therefore both particle and antiparticle excitations. For quantizing (4.104, 4.106) 



we need for this reason two copies of refiection-transmission algebras |63j. The first one 



is generated by {ai{k), a*{k) : A; G M} and §(/c) given by (3.49). We denote the second one 



by A\_ because its generators {bi{k), b*{k) : A; G M} obey the anti-commutation relations (2.1 
2.2) with the transpose scattering matrix §*(A;). Besides (2.3), one has therefore 



k{k) = j2b,{-k)s,,{k) , b:{k) = j2s,,{-k)b*{-k) 

i=i j=i 



(4.108) 



In what follows we use the convention according to which {ai{k), a*{k)} and {&i(/c), b*{k)} 



annihilate/create respectively antiparticles and particles. The solution of ( 4.104|4.106 ) in this 
basis is 

Mt,x,i)= — k(A;)e'*'=(*-") + 6*(-A;)e''=(*-")l , (4.109) 
""^ dk 



27r 



[ai{-k)i 



-\k{t+x) 



+ bKk)e 



ik{t+x) 



(4.110) 



In the Dirac junction the anti-unitary operator T of time reversal and the unitary operator 
C of charge conjugation act as follows: 



Tipi{t,x,i)T ^ = r]Ti'2i-t,x,i) , Ti)2{t,x,i)T ^ = r]Ti'i{-t,x,i) 



IVtI 



(4.111) 
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Ctpi(t,x,i)C ^ = -r]ci'l(t,x,i) , Cij2(t,x,i)C ^ = 7]ci)*2{t,x,i) 



\Vc\ 



1. (4.112) 



Like in the Schrodinger case the boundary condition (4.106) is invariant under time reversal 



only if U is symmetric (3.30). The condition for charge conjugation invariance is instead 

~ (4.113) 



U 



-U. 



The violation of (3.30) and/or (4.113) leads to the breakdown of the corresponding symmetry 



by means of the boundary condition (4.106). 



The electric current and energy-momentum tensor are 

3tit,x,i) = -e : ip*!/; : {t,x,i) , 
j^{t, X, i) = -e : ip*aip : (t, x, i) 



eu{t,x,t) = -:[r{dt^)-{dtr)^]:{t,x,t), 
e^t{t,x,i) = ^ : [ilj*a{dtip) - {dttp*)aip] : {t,x,i) 
where : • • • : denotes the normal product in A+ and A', and 



1 
-1 



(4.114) 
(4.115) 

(4.116) 
(4.117) 

(4.118) 



According to our convention for particles and antiparticles, the incoming asymptotic sub- 
algebra is generated by {ai{k), a*{k), bi{—k), b*{—k) : A; > 0}. The edge Hamiltonians 
and the asymptotic charge operators of particles and antiparticles are 



hi 



° dk 



k\[a*{-k)a,{-k) + b*{kMk)] 



dk 
2^ 



bUkMk) 



Qi 



° d^ 



a*{—k)ai{—k) 



(4.119) 



(4.120) 



respectively. We associate with (4.120) the chemical potentials /ij and /ij. Now, following the 



general strategy explained in section 2, we set 

n 

K = ^/3i{hi- HiQi - JliQi 



(4.121) 



i=l 



and define first the steady state ^sun on the sub-algebra by means of (2.9). Employing 



( 2.3[4.108 ), we extend after that the state ^i3,f_i,ji to the whole algebra generated by A+ and 
A\_. In this way one gets 



{a*{p)a,{k))p,,r, = 2vr{ \e{k)Mk)6ij + d{-k) U*z M-k) %! 6{k - p) 



1=1 



+ 



9{k)mv,, + e{-k)mf,{-k) 6{k + p) 



(4.122) 
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n 

1=1 



+ 



9{-k)v*Mk) + e{k)fj{-k)Vj, 6{k + p) 



where 



m 



m 



(4.123) 



(4.124) 



are the Dirac distributions for antiparticles and particles respectively. 

Let us discuss finally the behavior of ^i^^^^ji under charge conjugation. One easily verifies 
that ^/3,^,ji is invariant under charge conjugation, namely CVL^^^^ji = ^i3,^,ji, provided that both 
conditions 

fii = -Jli (4.125) 



and (4.113) hold. 



4.2 Transport properties 

We are ready at this point to derive the steady currents in the state ^i3,^,-ji- A computation, 
analogous to that performed in section 3, gives 



e 

2^ 



dk 
2^ 



1 + e^^-^^- 



_|_ Q-PjH 



f^{~k)-f,{k) 



(4.126) 



which satisfies Kirchhoff 's rule and vanishes at equilibrium exactly like the Schrodinger steady 
current. A new feature of the Dirac steady current is the presence of particle and antiparticle 
contributions, captured respectively by the numerator and the denominator in the fraction 



under the logarithm in (4.126) 



Concerning the dependence of (4.126) on the chemical potentials, some particular cases of 



are worth mentioning. We first observe that if charge conjugation is preserved, the steady 



current (4.126) vanishes because of (4.125). This is due to a cancellation between the particle 



and antiparticle contributions. If instead all /i, = the current (4.126) is temperature inde- 



pendent (in spite of the fact that the junction is in contact with heat reservoirs with different 
temperatures fij) and takes the simple form 



(4.127) 



In the case /ij = the steady current (4.126) coincides with that of the Schrodinger junction 



(3.50) with scale invariant boundary conditions. 
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Let us focus now on the temperature dependence of (4.126). At high temperatures one has 

(4.128) 



i=i 



whereas at zero temperature 



The conductance tensor corresponding to (4.126) is 

G.,(/3,/x) = |^(5.,-|U,,f)^ln 

and is not symmetric in general. 
For the energy flow one obtains 



7I(/3,/x,/i) = (^xt(i,a:,i))/3,M,M = {Sij - [U^f ) 



1 + e'^^^'^ 
I _|_ e-ft-Mj 



dk 



27T 



1 " 1 

^ E - ^ [Li2(-e/^^-^0 + Li2(-e- 



(4.129) 



(4.130) 



(4.131) 



One can easily verify that the charge and energy densities are obtained by the replacement 
(lUjjp — 6ij) I — )■ (|Ujjp + 6ij) in (4.126) and (4.131) respectively. Notice that these quantities 
are x-independent and therefore do not present Friedel oscillations. The reason is that both the 
dynamics (4.104) and the boundary conditions (4.106) are scale invariant for massless fermions. 



4.3 Noise 



For conciseness we report directly the zero-frequency noise power in terms of the matrix U 



appearing in the boundary condition (4.106) and the distributions (4.124). One has 



dk f 

— [6,,F,,{k) - \i],,\^Fu{k) - |%|2F,,(fc) + 

1 " — — 
+ - J2 U/,%U„,-U™[F;^(A;) + F^K^)] 



Lm=l 



(4.132) 



with 



F,,{k) = mil - f,{k)] + Mk)[i - f,{k)] . 



(4.133) 
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If all the temperatures are equal = (3), the /c-integration in the right hand side of (4.132) 
can be performed exactly and gives 
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(^..■-|u.,f) 



+ 



1 + e^/"' 1 + e"^^' 



l,m=l 



e/3Mi + e^^'" / 1 + e'^ 
In 



+ 



1 + el^^'^ 1 + e-(^^'3 
— m 



1 +e- 



+ 



. (4.134) 



For the purely thermal noise on gets therefore 



(4.135) 



which coincides precisely with the result (3.88) for the Schrodinger junction at criticality. 



4.4 Remarks 



In spite of the different dispersion relations, the general structure of the steady currents and the 
noise in the Schrodinger and Dirac cases are quite similar. A characteristic feature of the Dirac 
case is the possibility to introduce the independent chemical potentials /ij and Jli, associated 
with particles and antiparticles. This fact has elementary but important consequences. If 
/ij = —Jli the particle and antiparticle contributions cancel each other in the electric steady 



current (4.126), but sum up in the heat current (4.131) and in the zero frequency noise (4.135). 



Following the argument in section 3.5, the results about the Dirac junction have a straight- 
forward generalization to the case when ijj is minimally coupled to a static classical electromag- 
netic field generated by the potential (At(P), A(P)) in the ambient space. 



5 Outlook and conclusions 

In this paper we developed an algebraic method for constructing non-equilibrium steady states 
on star graphs. Our approach is microscopic and our construction generalizes that of a 
Gibbs state over the algebra of canonical (anti) commutation relations. The Schrodinger and 
Dirac equations have been investigated in this framework. We considered in detail the case 
in which the interaction, driving the system away from equilibrium, is localized in the vertex 
of the graph. It turns out that the non-equilibrium dynamics, generated by such interactions, 
is exactly solvable. In fact, the fi^^^-expectation values of various observables (currents and 
charge densities) can be computed exactly, without resorting to any kind of approximation. We 
have shown in particular, that the expectation value of the electric current in the Schrodinger 
case reproduces precisely the famous L-B formula. Once the formalism has been tested on the 
L-B steady current, we applied it for the computation of the charge and energy densities, the 



24 



energy flow and the noise power. The presence of Friedel oscillations has been detected. We 
demonstrated also that point-like interaction in the junction modifies the linear dependence 
of the thermal noise on the temperature ( Johnson- Nyquist formula). The formalism has been 
generalized in order to include the minimal coupling to an external time-independent electro- 
magnetic field as well. 

Summarizing, the star graph models proposed and analyzed in this paper represent relatively 
simple exactly solvable examples of quantum non-equilibrium systems in a steady state. For this 
reason they provide a nice laboratory for testing general ideas about non-equilibrium dynamics. 

Our results can be generalized in various directions. First of all, one can consider more 
complicated networks with several junctions and loops, which can be crossed by magnetic 
fluxes. The basic idea for treating this case is to replace in the above formalism the scattering 
matrix § with an effective one §eff, which takes into account all vertex interactions jMI-EHI and 



the presence of a magnetic fleld (see (3.101)). The derivation of the L-B steady current and 
the noise in this case is of particular physical interest and is currently under investigation |69j . 

Another possible generalization is the study of imperfect leads involving interactions with 
external potentials and/or self-interactions like those in the Luttinger liquid. More general 
boundary interactions, involving new vertex degrees of freedom of the type appearing in the 
resonant-level model [70], can be investigated in the above framework as well. We will discuss 
these issues elsewhere. 
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